Observation of anomalous diffusion and fractional self-similarity in one dimension 
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We experimentally study anomalous diffusion of ultra-cold atoms in a one dimensional polarization 
optical lattice. The atomic spatial distribution is recorded at different times and its dynamics and 
shape is analyzed. We find that the width of the cloud exhibits a power-law time dependence with 
an exponent that depends on the lattice depth. Moreover, the distribution exhibits fractional self- 
similarity with the same characteristic exponent. The self-similar shape of the distribution is found 
to be well-fitted by a Levy distribution, but with a characteristic exponent that differs from the 
temporal one. Numerical simulations suggest that this is due to long trapping times in the lattice 
and correlations between the atom's velocity and fiight duration. 
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Diffusion is a phenomenon encountered in almost ev- 
ery branch of physics. Its ubiquitousness stems from the 
central limit theorem, which states that a sum of ran- 
dom variables is distributed normally as the number of 
addends increases. It holds, however, only when the dis- 
tribution of the variables has a finite variance. When 
this assumption does not hold, i.e. for heavy-tailed dis- 
tributions with asymptotic power-law behavior with an 
exponent a -1-1 with < a < 2, the sum converges instead 
to a Levy distribution Lq,. A diffusion process which re- 
sults in such non-Gaussian spatial distribution is usually 
regarded as anomalous [1] . Heavy-tailed distributions are 
found in many fields, from animals foraging strategies [2], 
to the prices of commodities and stocks [3]. In physics 
they emerge in situations where there is no characteristic 
length scale, such as near phase transitions [4], in turbu- 
lent flow [5] , in quantum phase diffusion or when a sys- 
tem is out of thermal equilibrium ]7\. In fact, anomalous 
transport properties are intimately linked to non-linear 
chaotic dynamics which naturally appears in many phys- 
ical systems |H1 E] • 

A simple diffusion model we have in mind is that of 
particles in real space, each having a velocity which fluc- 
tuates in time due to interaction with a bath. After some 
time the particles' position is distributed as W{x, t). The 
characteristic width of the ensemble, e.g. the full width 
at half the maximum (FWHM) , usually scales as a power- 
law i^/". Anomalous diffusion can arise when the distri- 
bution of velocities has heavy-tails, and almost always 
results in a 7^ 2 [T]. The theoretical challenge is to 
connect the microscopic physics to the evolution of the 
distribution W{x,t). The experimental challenge is to 
measure these distributions in a well controlled and iso- 
lated environment. In this work we meet this challenge 
by measuring anomalous diffusion of laser cooled atoms 
in a polarization optical lattice [TUHlSj . In this system 
the steady state atomic velocity distribution was shown 
both theoretically and experimentally to follow a power- 
law, with an exponent that depends on the lattice depth 
[TTl fT2l [14] . It was also predicted that the real space dif- 
fusive motion of the atoms in such a lattice is anomalous 
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FIG. 1. a, The experimental setup as seen from above. The 
atoms are released from a small crossed dipole trap into a sec- 
ond far-off-resonance trap with a very elongated aspect ratio 
which confine their motion to a single direction (tube trap). 
This trap is superimposed with the two counter-propagating 
Sisyphus lattice beams which induce the diffusive motion. Af- 
ter a chosen time, the atomic density is recorded using ab- 
sorption imaging technique. Each image gives the column 
two dimensional density, and we integrate over the radial di- 
rections to obtain the axial atomic distribution, b. Atomic 
distributions after different diffusion times in a 4.8-Er deep 
lattice. Each measurement is repeated 20 times and averaged 
to improve the signal to noise ratio. 



for a wide range of lattice parameters [13] . The onset of 
anomalous transport characteristics was observed with a 
single trapped ion |15) . However, a measurement of the 
distribution W{x, t) and its anomalous dynamics was not 
reported to date. 

Here we report such a measurement in one dimension 
with an ensemble of ultra-cold ^^Rb atoms. By setting 
out with a very small atomic cloud and recording the 
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longitudinal density distribution after different waiting 
times, we are able to directly measure anomalous dynam- 
ics in the lattice. We find that the width of the distribu- 
tion exhibits a power law time dependence, from which 
a characteristic exponent can be extracted. The value 
of this exponent depends on the lattice depth. Further- 
more, we show that the density distribution at different 
times exhibits self-similarity with the same characteris- 
tic exponent. The self-similar shape of the distribution is 
found to be very well fitted by a Levy function. However, 
the characteristic exponent extracted from this fit is sig- 
nificantly smaller than the exponent extracted from the 
dynamics. We investigate this point using classical and 
quantum Monte-Carlo simulations and find that it origi- 
nates from correlations between the atom's velocity and 
the corresponding time it spends un-trapped in a partic- 
ular lattice site, combined with heavy-tailed distribution 
of the durations in which it is trapped by the lattice. 

The apparatus is depicted schematically in figure [l^. 
In each experiment, ~ 1.5T0® Rb atoms are prepared in 
a cigar shaped crossed dipole trap with an aspect ratio of 
1:3.9 and a radial oscillation frequency of 27r • 420Hz (for 
more details regarding the apparatus see Ref. [H]). The 
Sisyphus lattice is created by two counter-propagating 
lattice beams with identical wavelength and orthogo- 
nal linear polarizations [10] . These beams originate 
from an external cavity diode laser whose frequency is 
locked to an atomic transition and detuned — 66MHz rel- 
ative to the transition between states 5^Sif2,F = 2 and 
^'^Pz/2,F' = 3. The atomic cloud has its long axis aligned 
parallel to the lattice beams. The lattice depth is calcu- 
lated from the measured beam's intensity and waist of 
1.1mm, and using a saturation intensity of 3.6mW/cm? . 
Before switching on the lattice, the atomic cloud has a 
temperature of 12 and a maximum phase space den- 
sity of 1.7 • 10^'^. To counteract gravitation and improve 
the signal to noise ratio, the motion of the atoms is con- 
fined to the longitudinal direction by superimposing the 
lattice with a "tube" trap, namely a Gaussian beam with 
a waist of ^ 120fim such that the Rayleigh range is much 
larger than typical diffusion distances we record. Both 
the crossed dipole trap and tube trap originate from a 
single frequency Ytterbium fiber laser at a wavelength of 
1.06/im, but their frequencies are shifted relative to each 
other by more than 20MHz to prevent standing waves. 
The power of the tube trap beam is lOW. 

The lattice is turned on 1ms before the crossed dipole 
trap is turned off, during which the atoms equilibrate 
with the lattice. The average rate of photons scattered 
by each atom from the lattice is ~ 10'^s~^, much larger 
than the initial maximal average elastic collisions rate of 
190s~^. This is important since elastic collisions lead to 
an undesirable thermal equilibrium. We define t — a,s 
the time at which the crossed dipole trap is switched off 
and the atoms start diffusing in the lattice. The initial 
size of the cloud is ^ 200/zm, much smaller than the 
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FIG. 2. A log-log graph of the FWHM squared as a function 
of the time squared, for different lattice depths. To facilitate 
the comparison of the slopes, all graphs were shifted in the 
y-direction to cross at a single point. The linear behavior es- 
tablishes the power-law time scaling of the width of the atomic 
distribution. The FWHM is found by fitting the atomic dis- 
tributions to a general function from which we extract the 
width. In very shallow lattices (below ~ l.JE,.) we notice an 
excess density in the center of the distribution which is prob- 
ably due to corrugation in the tube trap potential. To avoid 
systematic errors due to this effect, we exclude in this range of 
lattice depths the central part of the distribution when fitting 
and extracting the FWHM. 



typical diffusion distances. We take a series of absorption 
images after different waiting times, an example of which 
in a A.SEr deep lattice is depicted in Figure [TJd. 

A convenient and useful theoretical framework which 
describes a broad range of anomalous diffusion processes 
is the fractional diffusion equation (FDE) JJj: 



dW{x,t) 
dt 



(1) 



where W{x,t) is the atomic distribution at time t, 
is the Weyl operator describing a fractional derivative in 
space and similarly D\~^ is the fractional time deriva- 
tive. Kp is a generalized diffusion constant having the 
dimensions of cm^^/s^ . For /i = 2 and f3 — 1 this equa- 
tion reduces to a normal diffusion equation, fi < 2 cor- 
responds to long spatial jumps (also referred to as Levy 
flights), whereas f3 < 1 corresponds to long dwelling times 
between jumps. The solution for the kernel Gp^{x,t) of 
this equation can be written in terms of Fox functions 
[18j . A general property of the kernel is its time and 
space scaling [T71 [THj : 

Gp^{x,t)=t-''/''Gp^{xt-^/'') , (2) 

where G^^ is the reduced kernel function. 

One conclusion that can be immediately drawn from 
Eq. [2] is that the typical width of the distribution should 
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FIG. 3. The difFusion exponent a as a function of the lattice 
depth. The exponent is extracted by three different methods; 
The (blue) squares are obtained by fitting the FWHM scal- 
ing to a power-law FWHM ~ t^^" (the errorbar represents a 
95% confidence level) . The (red) circles are obtained from the 
measure of the self-similar transformation. The (black) trian- 
gles are found by fitting the distributions taken after 30ms 
with a Levy La function. 



scale as t^^", with a dynamical difFusion exponent given 
by a = To test this, we plot in figure [2] on a log-log 
scale the FWHM extracted from the data as a function of 
time for different lattice depths. The curves are approxi- 
mately linear, indicating that indeed the width scales as a 
power-law in time. Furthermore, the slope of each curve 
is different, showing that the difFusion exponent depends 
on the lattice depth, as predicted jjjj il2j . We fit each of 
these curves with a line, and from its slope we extract the 
difFusion exponent. The results are plotted as squares in 
figure [3j and demonstrate that the whole range of frac- 
tional diffusion exponents is accessible in our experiment 
by changing the lattice beams power. Note that since in 
the radial direction there is no cooling, atoms are even- 
tually lost from the trap due to spontaneous emission. 
Owing to the large Rayleigh range of the tube trap, there 
is almost no mixing between the axial and radial velocity 
distributions and therefor and therefore the radial loss 
does not distort the axial spatial distribution. Neverthe- 
less, we use the data only as long as at least 10% of the 
atoms remain in the trap. Increasing this number to 30% 
changes only the 3 largest points in figure|3]by up to 25%. 

Another important conclusion from Eq. [2] is that the 
distribution should exhibit a self-similar scaling with re- 
spect to xt^^/°', regardless of its exact shape. An ex- 
ample of this property is shown in the inset of figure 
|4j When using the appropriate a, all experimental data 
taken at different times collapse to the same curve when 
the X and y axes are re-scaled according to Eq. [2j This 
property can be used to extract the dynamical exponent; 




FIG. 4. The measure of the self-similarity, m, as a func- 
tion of the diffusion exponent, a, for three dilTerent lattice 
depths. Each point is calculated from 13 distributions taken 
after 10 — 40ms of diffusion and normalized such that their 
integral is unity. The inset shows an example of the re-scaling 
transformation with a = 1.25 in a 4:.8Er deep lattice. 



we employ an Li-type measure of the self-similarity: 
JZ,\W^ix,a)-W{x,a)\dx 



m{a) 



E 



JZo ^(^' ")^^ 



(3) 



where W,{x,a) = ty^Wixt'/^^ti) are re-scaling of the 
series of distributions W{x, ti) measured at different 
times, and W{x, a) is their average. In figure|4]we depict 
this measure as a function of a for three lattice depths. 
For each of these curves there is a single minimum whose 
position is shifted in accordance with the lattice depth. 
We interpret the minimum as the most probable value for 
the diffusion exponent. This value is plotted as a function 
of the lattice depth in figure [Sj The diffusion exponents 
obtained by the self-similarity method and by fitting the 
FWHM to a power-law agree to within the uncertainty. 

Up to this point we have analyzed the temporal be- 
havior of the atomic distribution, and now we turn to 
study its shape. Motivated by the fact that for /3 = 1 
the solution for the kernel of Eq. [T] is the Levy stable 
law La [17l [T8] , we use the latter as a fitting function. 
In figure [5] we depict the shape exponent, a, extracted 
from these fits as a function of the diffusion time, for 3 
characteristic lattices. In all three lattices the shape ex- 
ponent converge to an asymptotic value on a timescale of 
10ms. The inset shows the fits after 30ms of diffusion. 
The Levy distributions fit the data very well with an av- 
erage r-square of 0.96 and an uncertainty in the fitted a 
at a 95% confidence level of ~ 3%. The shape exponent 
is depicted as triangles in figure |3] At very shallow lat- 
tices it approaches 2; this is expected since for any finite 
observation time there are lattices which are too weak 
for the atoms to reach equilibrium. In deeper lattices, 
on the other hand, we consistently find that the shape 
exponent is significantly smaller than the dynamical ex- 
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FIG. 5. The atomic distributions after different diffusion 
times are fitted with a Levy La function. The resulting ex- 
ponent a is plotted as a function of the diffusion time for 
three characteristic lattices. In the inset the distributions af- 
ter 30ms of diffusion and their corresponding fits are shown. 



and flight duration in order to sustain this result for in- 
creasingly deeper lattices. Both these factors were indeed 
found in the MCWF simulations [19] . 

To summarize, we have presented measurements and 
analysis of spatial anomalous diffusion of ultra-cold 
atoms in a ID polarization lattice. We find that a com- 
plete description of this process goes beyond the FDE and 
must include the effect of correlations between the mo- 
tion variables. Future extensions of our work include the 
study of anomalous diffusion in the presence of an exter- 
nal force and diffusion in very shallow lattices where there 
is hope to observe super-ballistic diffusion [5T]. Also of 
great interest is to measure the time-evolution of the 
velocity distribution which is also expected to exhibit 
anomalous diffusion [22l [23] . 

We thank Rami Pugatch for stimulating discussions. 
This work was partially supported by MIDAS, MIN- 
ERVA, ISF and DIP. 



ponents, and in particular smaller than 1. Similar results 
are obtained if the tail of the spatial distributions is fitted 
with a power-law instead of with a Levy distribution. 

In order to better understand why the shape expo- 
nent is different than the two dynamical exponents, we 
have written both classical and quantum simulations |19j . 
Based on their results we attribute the discrepancy be- 
tween the exponents to the combination of two factors: 
/? < 1 in the FDE and correlations between the parti- 
cle's velocity and flight duration. The quantum simula- 
tion we have performed is a Monte Carlo Wave-Functions 
(MCWF) of atoms in one dimensional polarization lattice 
with angular momentum Jg = 1/2 to Je = Jg + 1 = 3/2 
transitions [SH] . The atomic evolution in the light field is 
treated quantum mechanically and photon scattering is 
described by quantum jumps. In general, we find a good 
qualitative agreement between experiments and the sim- 
ulation despite differences in the level structure and lat- 
tice detuning which we introduce to simplify the numer- 
ics [H]. The temporal evolution of the simulated spatial 
distributions, as well as their shape and dependence on 
the lattice depth, exhibit the same properties that were 
discussed above. In particular, similarly to figure [3| we 
find that the spatial Levy exponent is consistently smaller 
than the dynamical exponent for deep lattices. The sim- 
ulations also establish a clear correlation between the ve- 
locity and flight-duration of the atoms. In other words, if 
an atom acquires a large momentum, it is more likely not 
to be trapped in a single lattice site for longer times. We 
have also ran classical Monte-Carlo simulations where the 
velocities, flight-durations and dwelling times in the lat- 
tice are all drawn from Levy distributions. Though this 
simulation shows that a shape exponent smaller than 1 
can be obtained by /3 < 1 for shallow lattices, it is nec- 
essary to include the correlations between the velocity 
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